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Preface 


The present thesis entitled 'Certain Investigations relating to blood flow 
in micro circulatory system' is a record of genuine research work carried out 

by me under the supervision and guidance of Dr. H.S. Shukla (Supervisor). 
Principal, Pt. J.N. College, Banda. 

The purpose of the thesis is to study some problems on steady and 
pulsatile laminar viscous flow of blood in micro circulatory system i.e. the flow 

of blood in those vessels whose diameters are less than 500 
(l^m = 10 m). 

The effects of peripheral layer, yield stress and other rheological 
parameters on blood flow have been analysed. We have studied the problems 
on dispersion process in blood flow and effects of mild stenosis on blood flow 
through narrow vessels. In the analysis it has been assumed that the radius 
of the vessel is constant and there is no slip velocity at the wall. 

This thesis consists of six chapters. Chapter I contains ctbe general 
introduction of the subject (viz. Circulatory system, rheology of blood, geometri- 
cal aspect of vessels etc.) and a brief acpount of the relevant literature. 

In Chapter 11, we have discussed steady shear flow of micropolar fluid 

through a channel of large cross-section. In the present chapter we have 

proposed to analyse the fully developed steady shear flow of micro polar fluid 

between two parallel plates. The polarity arising due to the presence of buoyant 
corpuscles in the fluid leads to give a theoretical insight of the blood flow 
through large arteries. The microrotional field depends on the length ratio 
parameter. Microrotation effects and microrotational inertia are considered in 


this case. Exact steady solution for velocity and microrotation subject to the 
strong and weak limit of influence of surface on microrotation are obtained. 

The parameters characterising the ratio of vortex viscosity to sheer viscosity 
and length ratio which do not appear in a classical Newtonian fluid are 
introduced and its effects on velocity and microrotation are discussed. The 

results are discu.ssed with the help of Tables. 

In Chapter III we have discussed steady flow of blood in Narrow tube 
through r.Hcropola^uid model with slip at the wall. In the present chapter >5 

we have studied the poiseville flow of micropolar fluid as a model of blood 
flow in 40 ^m diameter tube and for 40% RBC concentration. 

At the boundary, we have introduced a slip in the axial velocity and 
a partial rotation depending upon the wall effect paramters and fluid vorticity. 

From the analysis we have observed that the introduction of slip reduces the 
apparent viscosity of the suspension and increases its axial velocity whereas 
it has no effect on particle rotation. Effect of S (0 ^ S ^ 1; a parameter 
depending on S and concentration) on different flow parameters are seen and 
by comparing the results is in the experimental values we have tried to estimate 
the value of S. Results have been discussed with the help of tables. 

In Chapter IV, a mathematical analysis for the dispersion of soluble 
matters in blood flow has been carried out. In present problem we have prepared 
the flow model with casson fluid in the core region surrounded by a Newtonian 
plasma layer near the wall. Taylor's limiting condition and Pick s Law of diffusion 
are used for finding the solution of problem. The effective dispersion co-efficient 
with which the solute disperses across a plane moving with meatx speed of the 
medium is found to be decreased with respect to the yield stress and molecular 
diffusion co efficient where as a reciprocal behaviour is observed with respect 
to the viscosity of the casson fluid. The results arc discussed with the help 
of tables. 

In Chapter j3V, we have discussed the effects of wall layer thickness on 


blood rheology. From several studies it has been established that for blood 
flowing in narrow tubes (500 ^m to 50 ^m) a cell free plasma zone exists near 
the wall and the extent of this depends on many rheological parameters. From 
the analysis presented here, we see that apparent viscosity of blood decreases 
as the thickness of plasma layer increases. 

In Chapter VI, a theoretical model of steady blood flow in narrow vessel 
have been considered. In present analysis we have assumed that fluid in core 
region satisfy casson equation and in marginal layer region satisfy Newtonian 
equation. Apparent viscosity of blood has been determined as a function of yield 
stress, vessel diameter and peripheral layer (PPL) thickness. 
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GENERAL INTRODUCTION 


1.1 INTRODUCTION 

Bioengineering is the field which has been 
developed recently by cross fertilization of ' the 

■ ® 

engineering and biolpgical sciences, so that both are 
going to be utilized more iffectively for the benefit of 
mankind. 

Biomecihanics is an important branch of bioengi- 
neering. This is the study of the mechanical aspects of 
the* structure and function of biological systems. 

' Biqfluid mechanics an important branch of biomechanics# 
is the study wherein'.the principles of fluid mechanics\are 
ap^T ied' to understand the problems in biology. Fluid 
mechanics is essentially the science concerning properties 
and motion of fluids. In the presei^t ’" . we V 

shall concentrate only on the application of fluid t 

mccha||ics to study tho hlotxi rheology and blood flow in 

ii ' ^ ^ " 

human circulation . - 

'Blood', the red fluid of the blood vessels 
with which we are all familiar, is the transport] medium 

of the body. It is the medium by which all | living 

tissues are related to their external environment^ 

i.e. from the outside world. Blood is pumped to and 
from between tissues and elementary canal to take up 

oxygen and nourishment respectively, and it helps 




I 
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in returning the t^roducts of oxidation and metabolism 
in the tissues to the outside world by. the lungs, the 
skin and kidneys. It is also the medium by which growth 
and repair substances are transported to the tissues 
and by which the various controlling glands of the body 
can distribute their cliemical ^ uiessenyers . 

Blood is some what viscid fluid. in man and 

in all other vertebrate animals with the cxcfjption of 
two (-the amphioxus and Icptocoplialus ) , it is i od in colour. 
It consists of a continuous yellowish aqueous phase, 

the plasma, in which formed elements are suspended - 

The formed elements consist of red blood cells (erythro- 
cytes), white blood cells (leukocytes) and platelets 
(thrombocytes). Plasma is made up of water (92%)‘ 

and contains traces of l^iorganic and organic salts. 
The inorganic and small organic molecules contribute 
little to plasma viscosity which is primtiiily dependent 


on the plasma proteins 


If blood is allowed to clot 


and the solid material be removed, the remaining fluid 
is called serum. Tliis lias essentially the same 
as plasma except that fibrinogen and t-'rae clotting 

factors have been removed. About 7%, by weight, of plasma 
is protein mainely albumin, globulin and flbiinogen, 
which have molecular weight ranging from 44,U0() Lod0,00000. 
About half of the protoir, mass is albumin. The 5igni.ficance 











of plasma procein is its multi form natures mentioned 

T 

below :- • ' 

i) Plasma proteins are responsible for osmotic 
pressure, the level of which is impori^ant for 
regulating water exchange between blood and tissues. 

ii) .They possess buffer properties ami maintain the 
lac id “base eguilibri'.im <,)C blood . 

iii) 'They produce a definite viscosity of plasma which 

is important in maintaining blood pressure. 

>iv) • They promote stabilization of the blood by providing 
conditions that prevent sedimentation of 

erythrocytes.", ' 

v) They play an important role in coagulation. 

vi) Plasma proteins are important factors in immunity. 

If blood, to which an anticoagulant has been 

* 

added, is poured into a test tube and centrifuged, the ^ 
corpuscles tend to sottie at the bottom and the blood 
is divided into two layers, namely, a red layer lower 
consisting of the formed elements and a transparent colour- 
less or slightly yellowish upper layer consisting of 
plasma. The leukocytes form of 'thin white film between 
the erythrocytes and tlie plasma since th^r specific 
gravity is less than that of the erythrocytes. 

I 

The red cells surface carrier a negative 

charge, but in stationary blood, cells interact with 
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each other to form aggregates. The aggregates commonly. 
«? 

consist of 6 to 10 cells stacked face to face and such 
a cluster of cells is called a rouleaux. Secondary aggrega- 
tion of rouleaux also occurs building up a coaplex three 

dimensional structure. When blood is sheared, fc-hose rouleaux 

breakup and at sufficiently high shear rates the cell 

gjjists as an individual. Blood flows in the form of 

different laminar containing different types of cells. 
Considering these factors, the rheological ^-'roporties 
of blood might be expected to be rather complex. 
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The red cells out number the white cells and 
platelets by approximately 5000 : 1 and 12 : 1 respectively. 
In noriiial blood, cellular components other tjjati red cells* 


I 




compose onTy about 3% of the total cell* volume ^nd usually 


have only a small effect on blood viscosity. A suspension 


of 

Vdlito cell 




Hi 

DHol 

Further, 

V- 1 

. 

' f 

blood 

viscosity 

X' ‘ 

^ ■ 

in 

leukocytes 



concentration (Dintenfass [38] ) 


r--' 


Tne red blood cell? are about 911 of, the total 

0 

cell volume in the blood. Con'soguen Ll y , removal 

of white blood cells and platelets does not measurably 
modify the experimently determined flow properties of 
these suspensions. Tlie concentration of mo^tr red blood 




volume fraction occupied by the red blood cells, called 
hematocrit. The hematocrit is normally about 42-45% 
of volume. The erythrocytes consist of a thin flexible 
unstretchable membrane with an interior filled with saluratoj 
hemoglobin solution (viscosity 6.0 centipoise). The 

membrane is higl^ly doforiuablo lJ the cliaiujo in the iniLfacc 
area is small, but becomes very much stiffer if the 
deformation produces a largo change in tho aioa t.he 

cell membrane. The red blood cells have the shape of 

a bictoncave disc which can deform, however, into a bullet 
shaped entity during passage t^irough small capillaries* 
The cell has a large surface area relative to its volume. 


Blood remains fluid in the blood vessels t'l t. oughout 
its life but rapidly .bpcomes solid when shed. These 
two qualities are essential for the preservation or life. 
The maintenance of fluidity is necessary lor ’Llio cii:„uiation 
of blood, whilst the solidification of the red blood 
provides an indispan_sable defence against excessive bleeding 
from wounds. The coagulation of the blood is due to 
formation of a ^elly by the deposition of protein matjerial 
called 'fibrin' and it is tlie formation of this body 
that the fundamental change in blood clotting occures. 
To prevent tliis from occuring, various anticoogulfehlTS ; 
are added to the blood when it is drawn from tlie' animal. 









The effects of an anticoagulent on rheological properties 
of blood are difficult to as^j^ertain because of the di.fficulty 
in making reliable rheological measurements in t^e short 
time available between the blood drawing and the blood 
clotting. Tho effect of apti coagu 1 en t-H on blood viscosity 
appears to be small (Cliartn and Kurland [20) ; Cokelet 

et al p4j ) . j ■ 

. f 

/ * 

Human Circulatory System : 

The circulatory system of blood consists of complex 
functions (Respiratory, Nutritive exetretory , ‘rotectrve 
.and regulatory) in the human body. The provides 

the energy for the circulation has four ^compar tmeiit {i*e* 
left and right atrium and left and right vorj. rides ). 
interconnected to each other by one way vaiveS* Blood 
comming from body tissues enters the right atriujTn through 
the venui cavul contraction of right atrium trCS'ii^eS blood 
past the Tricuspid valve into the right ventricle. From 
this point, there begin two subdividsi^pns , viz. the pulmonai 
and the systemic circulations. As the term shows, the 
former services the . lung and^ latter the various systems 
of the body. The blood vessel has internal diameter 
in the range of 2.5cm in. aorta to about 4 microns in , 
capillary (Keela and Neil [ 7/7 • 

Blood is pumped from the heart into the aorta 
from where it goes to the circulatory system tdnsisting 
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m 
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1800 terminal branches, 4,00,00,000 arter 

1,20,00,00,000 capillaries, 8,00,00,000 venules, 1,800 

terminal veins, 600 main venous branches, 40 large 'venous 

branches, 40 large veins, and then returns to the heart 

« 

thro'ugh the ven>ule caval 

It is estimated that the heart muscles itself 
consume about 10% of the energy reguired isustaiu life* Only 


some . energy goes into the mechanical work of pumping 
blood. The heart beats about 70 times in onp minute 


in an average person at rest 


Blood flow . in smaller 


blood vessel approaches .steady flow condition other wis< 


vessel 
iie. T 


pulsatile. The pressure in the aorta rises rapidly to 

its maximum (systolic) value of about 120 mm ^Kg . The 

requirement of the circulation is the supply of _ oxygen 
required for ‘metabolic pressures. In an average man 
at rest, the O requires 200 ml/iuin. Under physicul stresse^s, 
the need may rise to 5 lit. /min. Blood has a capacity 
of 200 ml of oxygon per litre.; The efficiency oi transfer 
of ©2 from blood to the cells is such that for an average 
man it is necessary for the heart, to circulate 5 to 6 

litres of blood per minute or about the 5 times tiie metabolic 
requirement of oxygen. The 'metabolic j'equirements ' 

are the energy requirements Cor biological Cunct i c;n:-: . 

In fluid mechanics, the circulatory system is 

classified into two parts, viz. macrocircu la tion and 



microcirculation 


(a) Macroci rculation 


The flow of blood in vessels of diameter 
greater than 500 (such flow occures in aorta) 

is called macro circulation. The flow is characterized 
by high Reynolds number# defined as the raj.io of 


inertial to viscous forces 


Turbulent flow can 


occur in blood vessels for tl»e Reynolds number greater 


than 2300 


The governing equations for analysing 


of such flow conditions should 


includs* inertial 


effects, effect of curvature of blood vessel tjulsatile 
flow and the di s tens ibi 1 i ty of the voi\wcl wail. 
Blood can be taken as homogeneous and continuum 
fluid When vessel diameter is largo in c' xiipa r i s ion 


to t.ie dimensions of rod bipod coils 


I'u Lsatiic 


Clow effects and iJirciiMurc wave i imis have 

been studied in^ macrocirculation by a n.ii\ber of 
research workers, SWJilak and Tuyioi [iSl] . 

Tl\e flow disturbances at bifurcations, bcncis etc. 
and their effects on pathological states aid studied 
by l‘j^tel et al 

Microcirculation : 


When the diameter of blood voss(gi is . loss 


than 5 0 0 ,.Uin ( i . c . 


J. aLteiiolcs, capillaries, venules), 


circulation is calli'd 


m i c r o c i r c u 1 a tt'a vi' 
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is responsible for 80% pressure drop in circulatory 
sysjtem. In capillary bed the transfer of nutrients 

II 

to 'and removal of wastes form the living cells of 
the body is a part of microcirculation. The flow 
is characterized by very low Reynolds number's (Re 

I 

<1). But wo can not neglect the size of tti'e red 
blood cell, compared to blood vessels size, it be- 
comes necesfeary to treat the flow as two phase flow 
(:.on homogeneous). This two phase natur'"' of flow 
is responsible for anomalous behaviour in blood 
flow through narrow tubes. The diameters of the 
capillaries ranges are approximately 3.5 ^m to 

10 ^m. The problem is generally analysetl by either 
higli speed photography or by model studies. (Branemark 
[10], Gross and Aroesty 


Here we are interested in microci cculatory 
system . Table 1.1 gives ba.s.ic lnt-onii4bickii) related 
to the properties of the vessels through which the 
circulating blood passesr When blood fio.v through 
capillaries of internal diameter equal to tliat of 
the red cell, the flow pattern is entirely different 
from that in the large blood vessels. The red ' blood 
squeezes through the capillaries one after the other.- 
with plasma trapped in between them. Hegen-poiseuille 
law, axial accumulation of cells and the parabolic 
velocity distribution have not relevance tfc' capillary 









flow. The flow pattern, dovdloped when a colls 

follows each other throuyh such a narrow tube, fias been 

tortucd as bolus flow. Between the rear end ot one cell 

and the leading end of the next, circulation patterns 

are set up in the axial core of plasma. The efcect of 

this flow is that abolus of suspending liquid between 

two cells is carried down in the tube at the velocity 

« 

of the cells. 


1.3 


KIIKOLOGY OF BLOOD: 


Blood behave like a time dependent non-Newtorian 
fluid and the basic rheolog leal property of ‘uiood is 
its viscisity. The viscosity of whole blood is about 


4.0 centi.pose and of plasma is about 1.2 centipose W-lerril 
[^9^ 37 °C. The specific gravity of whole blood 

is 1.05 to 1.06 (Cokelet , [34] ) and that of plasma is 
about 1.03. Thus the red blood cells tend to sediment 

4 , 

slowly in plasma. /The increased viscosity relative 
to water is produced- mainvTy by the presence of plasma 
proteins. Both, the molecular shape ai)ci concentration 
of the protein are important. Fibrinogens which .has 
an elongated molecule has a marked influence; although 
forming less than 5% of the total plasma protein forming 
less than 5% of the total plasma proteins- it is res- 
ponsible for about 20% of the plasma viscosity elevations. 


1.3(a) The Vlscosityi of B lood: 

/ 

It JLS an important factor ^ determing the local 
pressure variation tiirougli cardio-vascular systc.ii , which 


X 






f'i - 

• /"V ; 


‘ 1 ^ 


r 




V, ■ ■( Y \ 




in turn influences the localf low rates through eacft^' section 

of the vascular network. The clinical imjjortance of i blood 

viscosity as a parameter lies in its scnsitivi Ly ,^to small . 
variations in composition. One can often diagnose, patholog- 
ical .states by determining a change ih'blbod viscosity . 

Thereareseveralrhcpiogical parameters I c.g. plasma, 

I blood cells, hematocrit etc.) which effect the blood 

* viscosity. The viscosity^ of plasma increases with ^ its 

f ' " ' ' 

protein concentration. But some proteins have different 

1 influences on plasma viscosity depending on their shape and 
siae. The influence of Fibrinogen on plasma viscosity can 
be' seen in the difference between plasma and serum viscosit- 
ies. Serum usually has a viscosity which is 20‘< less t^an 
"• . ■ ' ! ' ' ' *■ ' ' ' 

i that of plasma. The influence of globulins on viscosity is 

^ I ' '■ , ■ “■ - ^ n 

illustrated the disease macroglobuli^hemia (Sener [135J * i 

I ^ Albumin, th6 smallest plasma protein molecule ITable 1.2J 
is present in the largest concentration* Changes in albumin 
>; structure have tlio least effecL .of the protein oii plasma 

j; " viscosities, but the substance makes an important pontributt 

I ■ 

ion tcL plasma viscosity through its high concentration 
in plasmn* Many relationships have been suggested to expres" •>k 
i blood viscosity as a function of cell concentration, plasma 

viscosity and sliear rate* Wlion tciuporaturc is inci'^ascd# th© 
viscosity of blood and plasma is fallen. Measuremenis should 
be made at constant temperature ( 37 °c ) . : 

i Platelets and white cells/ in general, have little 

inf luences on blood viscos.Lty , because Lhcy presen t 
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in very minute quantity as compared to .red cells. j The 
deformation of red blood cells allow . the blood to remain 
fluid upto hematocrit of 98% j rigid cell without deformation * 
will cease to flow at a cell concentration of about 60%. The 
. fa^l in viscosity with increasing shear rate is also a 
f amif estation of the .deformatiop of the erythrocytes (Murata 
The hematocrit also ini'lLiicnco tlic do! onna L ion of 

the erythrocytes i.e. raising of the cell concentration 

/ 

produces an increase in cell deformation (Goldsmith [BoJ 
and therefore fall in viscosity of blood occures. 

In stationary blood, rouleaux is formed which Intracts to 
produce larger aggregates . At low flow rates the presence 
of these red cell structures strongly influences the 


viscosity of the blood. 

The size 

of the 

roule.iux 

and 

aggregates progressively 

decreases 

as the 

she- a r 

rate 


increases and this produces the typical shear thinning 

seen in blood at low shear rate. In normal blood the disa- 

ggregatioli is probably coiiipleto at shear rate 5*0 scc”'^ 

approximately (chien [33']). The viscosity of blood increases 
v/hen the shear rate falls, but it is uncertain as to what 
happens when the shear rate actually falls to zero. 1 The 
build up of a 3-d of increasing aggregates sugu ^ats that 

blood may, show a yield stress. 

l‘3{b) Blood as a Non Newtonian Fluid ; 




Blood, from fluid mecluuiics point of view, can be 
considered as a neutrally buoyant suspension of erythrocytes 
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in a Newtonian lisiuicl called plasma in small vcpsels 

the dimensions of red blood cells are no longer negligible 

I 

as compared to the tube size, hence blood can not be considr: 

ered as -a homogeneous fluid in these tubes of diameter of 

less than 500pm. The flow of blood througli this type of 

tube is of physiological and clinical importance. Due to 

its complex and anamalous behaviour it is very difficult to 

analyse it. The anamalies include Fahraeus-Lindqvist 

effect, segre silverberg effect, slip at the wall etc. This 

is specific nature of blood that exhibits a yield stress. 

Hence if applied shear stress is below a critical value, the 
/ 

response will be elastic and cn removal of stress uhe shape 
of blood film is unaltered. Otherwise flow takes place and 
blood behaves as a non-Newtonian fluid. At higher values of 
shear rates the blood shows Newtonian characteristics with 
a limiting viscosity. 


Some experiments by Bugliarello et al [16_] and 
Cokelet [134*J show that blood flow may deviate sigu.i f icantly 
from the Newtonian behaviour. Also it is found- that suspended 
blood cells are responsible for the non-Newtonian nature of 
blood rheology through . such mechanics as erythrocyte 
aggregation (coll protein .interaction) and erythrocyte 
deformation. (cell-cell and cell plasma interaction) (Chien 

#■ 

et al yij ) . Anamalous behaviour of blood show that ^at low 
shear rates the blood exhibits yield stress and bel'-ives like 



14 


suyyested H€.vSchc»l has al<u> ifceoi used to 

describe the Shea'r rate depencieuce oL blood. Charm and 
Karland | , sh&wcd a linutecj Shc=dr.- rate eha.Kje-. 

lida and Murato have studied the pulsatile flow 

of blood Lhroiujl, narrow riyid circular tube u.sincj HerschoJ- 
Bulklcy^fiuid model blood. I^uyliarello and Sev.lU|lU| have 
assumed J Casson fluid model -of blood in the study ar.d 
pulsatile blood flows in fine ylass. tube. Arimah et al ,3, 
have studied the similar problem thrdugh small riyid circular 
tube usincj micro-continuum api^roch and empto/)ry the 
condition of cell rotation at the solid boundary. Prahlad 

and Schultz yave two fluid models for hlood flow 

thiouyh liiuciii diviinotor Lube. 


Blood rheology can be hclplul in diagohosls of some 
blood abnormalities (Herrill ^ detreane in ' 

'mentation rate, increase in normal blood vjficotify and 
yield stress would indicate I'oly thymia . Anaemia ,s‘''usually 
indicated by lowcri.uj oC yield stress and in;o„uencemor,t o£ 
more Newtonian behaviour. L'levation o£ yield stios.s and extev, 
of non-Newtonian regime beyond 200 Sec'^ would indicate 
Hyper f ribrinoyonimia . 


1.4 


SUhVIlY OK l.rnJUAiiJHlt ; 


The German form of thd word (Rheoloyic) atxl the de- 
scription of a small vi;>comei:er as a ov" arc 

found in literature since early days; but, in its .modern 
senue, the term Rhcoloyy wd.s coined by Proi . Binyi.am 




■f' 


and korinally adot>l;Gd and defined at the foundaHon meeting 

of the American Society of Rheoloyy in December . 1929 

in Washinyton, D.C. as -The Science of the deformation 
and flow of matter". 


However, the sub;^ect of hydrodynamics and 
aerodynamics arc not included in rheologyi, Many authors 

(t^'liovre ot al Nubber etc.) have made studies 

on blood rljeolocjy . 


^ Records show that the Greeh and Romans were i^ttifre 

familiar with the anatomy of the heart and major i)Lofcd 
vessels. Roman writinys indicate that tho values of 

the heart were known to permit flow in only one iireotron 
The knowledge concerning the mechanics of circulation 
remained at the ievel of the ancient Greeks and Romans 
untll{_the renaissance when studies in anatomy led the 
way to the modern scientific approach. 

The most remarkable observer . Leonardo da Vinci 
cewbihed the study of structure with the 



i Oil the •Ccii.'tii O- Vascu J a r 

SySi-ems, Huch 

Of 

II 

his 

thinking 

reflected the classical 

schemes, iVn i ud i n g 

the 

idea tliat 

tile blood lias an ebb and 

flow rriotioft. 

He 

also 

remarked 

on the thickeniny and 

the hardening 

of 

the 

arterial 

walls, which we now call 

atherosclerosis 

« 

The 

modern concept of the circulation of blood was given 

by Willj.am Harvey in 1958. 
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The first mathematical paper of bloc<i flow was 

9 

given by Leonhard Euler in 1862. He deveiopal the one 

dimensional equations for inviscid flow of an inc ampressible 

! 

fluid in an elastic tube. His equations included both 
the conservation of mass aiid equations of motion. He 

also postulated a non-linear law relating the pressure 

» 

at any point inside blood vessel to its cross sectional 
area. Euler recognized that the heart could be thought 

of as a positive displacement pump and he gBve a fairly 
complete statement of the governing equations of blood 
■flow in the arteries. Hut lie could not solved the problem 

a 

posed by him. Recently, . many workers have studied the 
*flow property of blood (Bloor fS] , Brunn Copley [35], 

Goldsmiufi i j't. Liepsch Sugihara ! ■ i l, i, ' et v’ . 

Blood flow in capillaries is of great interest 

to i-ihys i o J og i ats involved in micro-vascu 1 u r research. 
The first direct observation of capillary vessels was 

reported by Malpight [ j | in 1661. In i-he beginning 

of 19th century, tlic principal anatomical features such 
as flow velocities, pressure and viscous loss phenomena 
were known quantitatively as \;oll. During igth century, 
general knowledge of the circulatory system Jjecome much 
more detailed and quantitative. 

^Ihe pressure drop relation in micro- c i rculati on 
was obtained by poiseuille in 1840. After t.his there 
appears to have been no further advance unt.ill recent 
times i'ii. -t-hooreticnl mechanics of viscous effects in 
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blood flow. Tlio most romai;kable discovery wa.s of Fahraeus 

# 

and LindciuisL 1'^ 'I ] / sliovv^lncj that the atJtjarent viscosity 

of bipod decreases iis tubi' diameter decreases from SOOjjm. 

f 

This also has been confi rmed by other investiqators I Coklet 
[]34] / Dintenfass 1*39] ) . They observed experimentally 

that this effect is seen to occur with capillary radii 
of as little as 5;ini, but this. Is reversed and the; apparent 
blood viscosity increases with a further decrease in 
the ca pi 1 a r ry radius. Recently, apparent: vi.sco.sity and 
otper properties wore studied by Copley [35]. A variety 
of investigations were proposed by authors such as study 
of capillary wall, the elaboration of Pick's law of 
diffusion and its application to. physiology etc. In 
biological science, tlio study of diffusivity of neutrients, 
me t. a l)e> 1 i i’ luuduelfi, drugs and other solutes is of at tno.s? 
importance!. Specially, many life giving materials mixed 
in the blood roach different parts of the body by the 
process of diffusion. Taylor studied' dispersion 

process in Newtonian fluid flows through a c Lrcular tube. 

<f Patel and Sirs (^hqJ and Rudraiah et al QxjJ made significant 
contributions in this direption. 

1.5 Geomotrical Aspect of Vessels : 

A striking characterxStic of the circulatory system 
is its geometrical complexity. Blood must flow through 
many bends, bifurcations, stenosis and tapering during 
its journey. While the complex anatomy of the major 
vessels has been known since ancient time, most of the 
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liiOO i'eif.s i n> wtjdi* rnta nd i luj the media ii i Ci3 of the c i rcu lattOt") 
is based on the inves tiyations of flow in u n j f csf jrrf 

'Cubes. It is only within the last two deciides that the;. 
effects of geometric transitions on blood flow have 
begun to be explained. 


1 5 ^ ^ S tonos i s : 


e 

A ijt'onotilH Ciluijefj tlio narrowing o-f- the blood 
vcssles due to the development of abnoroai tissues 
and gives way to serious circulatory disorder by reducing 
and oci^eluding tlic blood supply. Hemodynamic ch.a rac teres t ics 
may Jx) changed this undersirable growth which could be 

injurious to normal health. in arterial system of 

mammals, it isguite common to find narrowing?, .-some .of 
which avo at least, approx iraa tely , ax i symmetr i c; . These 
narrowing may bo caused by i n t r a voscu 1 a r or ,.0. 

implingcmcnt of i ioumehtn^ Or ^iguares on i hu vessel wad 


(Roach Rodbcird fi': 


[;ri-oj 


In recent years many workers have investigatec 
tiui flow dun.-acLoros L I c; ul blood thrcjugh ar'tery Hm 

prescncc of mild stenosis (Abdallah and Hwang f , 
Chakrawarty Malirola-a and Jay rarnan r g,, - c-'C.i. 

Young rf{>Mjgavo a theoretical analysis of th.e 
effect of time dependent stenosis on flow char acterestics 
of blood. He concluded that resistan.- to flow 
(impedonce) and the wall shear stress increase with increase 
' in the stenosis si/,e. Various authors |.avo scudiecl 
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Pulsatile efl.ect ol; blood in mild stenosis. Due to i.»resence 


of stenosis in the lumen, higher resistance to flow 
is caused and this effect becomes more prominant as the 
size of the stonos.s increases, also shear stress increases 
with respect to stenosis height ( Srivas tavi- 
Young CIS 41 . 


The geometry of the constriction choosen was 


defined by a' Guassian normal distribution c*jrve-. 


>evora . 


investigators have considered the numericai Lreatmout'v 
of the symmetric constrictions in rigid tw) dimensional 
conduits. Mac Donald C-O ^ Morgan and Young CSS") 

Rand et al- ['•I'T'T considered the non-Newtonian behaviour 
in their studies. Recently, Sinha and Singh {|5o3 studied 
the , effect of stenosis on blood flow through the coup^le 
stress fluid model. 


1.5(b) 'IT 


The mammalian arterial consists of the n^bin aortic 
tube which continues into the two iliac and femoral 
arLei |i.?s and a number of ujde branches. Transformations 
and propagation of pressure and flow waves depend on 
the distribution of the characterestic impedence. The 
two contributing factors are the decrease of cr'jus sectional 
area and the peripheral increase of the wall stiffness 


(elastic tapering) 


blood vessel.? have the 


clia iiic; 1 o iL'.s L i i: of "itcipor i;y 1 tndorS ratlier than straight: 


■ cylinders. This idea was given _by: Block i_8j • He gave 
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the idea that the vessels which carry blood tov/ards the tissues 
^ should be considerd as long, slowly tapering cones rather than 

cylinders. He obtained only pressure flow relation. The 
relations for shear stress and their variations with suspensiosn 
concentration and tapered angle, which, are also important, have 
been discussed in this thesis. 






TA B LI;: 1.1 
BASIC i.NKIBMATION 



Av . diam. urn 

Rate of Shear 
at wall Sec”'^ 

Re. number 

Large arteries 

10,000 

1,500 

5,000 

Terminal arteries 

1,600 

2,300 

200. 

Arterioles 

40 

5,000 

0.25 

Capillaries 
\ * '« 

16 

2,000 

0.02 

Venules 

60 

2,100 

0.25 

Main Veins 

4 , 000 

1,600 

800 

Lar<jc Vc ins 

20 , 000 

1,500 

20 


. » TABLE 1.2 

PROTEIN CONCENTRATION IN HUMAN PLASMA 


Cone . 

(g/1) 


Molecular- 
Weight (iftl) 


Mol- Sizes 
(uia) 


Protein 
A 1 bum in 
-•iobulin 
fibrinogen 


69,000 

1.50.000 

3.40.000 


15x4 
24 A 4.5 
70 X 4 


B 


P 

i 

1-*^ 

I 

1 


■ 


I 

iB''' 

I' 

1 

I 





I 
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eHAP<lER S7Sr 

Steady shear flow of micropolar fluid through a 
channel of Carge Cross Section 

introduction : 

The classical Navier - stobes theory does not adequately describe the flow 
properties of physiological fluids, e.g. blood, synovial fluid, etc. Mathematical 
model for the description of such fluids, which exihibit certain microscopic 
effects arising from the local structure and micromotion of the fluid element, 
is categorically introduced by A.C. Eringen (1964, 66). A subclass of these fluids 
is the micropolar fluids which have the microrotational effects and microrotational 
inertia. Ariman et al (19 74) reported that the theoretical results of polar 
fluid have good agreement with the experimental data of steady and pulsatile 
blood flow given by Bugliarello and sevilla (1970). 

In the present study we have proposed to analyse the fully developed 
steady shear flow of micropolar fluid betwcet two parallel plates. The polarity 
arising due to the presence of buoyant corpuscles in the fluid leads to give 
a theoretical insight of the blood flow through large arteries. The microrotational 
field depends on the length ratio parameter. Microrotational effects and 
microrotational inertia are considered in this approach. Exact steady solution 
for velocity and microrotational subject to the strong and weak limit of 
influence of surface on microrotation arc obtained. The parameters character- 
ii^ing the ratio of vortex viscosity to shear viscosity and length ratio which 
do not appear in a classical Newtonian fluid arc introduced and its effects on 
velocity and microrotation are discussed. 
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JAathematical analysis : 

Consider the constitutive equations for micropolar fluid in the form 

Tjj = (-p + ^ d|^j.) 5ij + 2p djj + 2 Hi ejjj. (Wj^ - aQ (1) 

~ ”p,p ^ij + P tT'i j + y CTj j (2) 

Where -Cjj is stress tensor 
Mj| is couple tensor 
Ejji^ is alternating tensor 
a^. is microrotation 
W|^ is vorticity 
p is thermodynamic pressure 

dij ^ y. (Vij + Vj_i) 

Wk = r.kij Vj i 

(I and p.’ are the viscosity coefficient of classical fluids and a, p, y 

are viscosity coefficients of micropolar fluids. These coefficient are subject to 
the restruction 

P] > 0, p > 0, y 0 

3p' + 2p > 0, 3 a + P + y ?! 0, -y <; p 5 y 

Incompressible flow of micropolar fluid between two fixed parallel infinite plates 

at distance h apart is considered. The flow throughout the region is maintained 

by a constant pressure gradient d^ . The lower plate with origin is assumed 

dz 

to occupy the plane y=0 and the upper plate is at y = h. At sufficiently large 
distance from the entrance region, the flow is supported to be fully developed. 
In the steady state the velocity v and microrotation c at any point have the 
component (0, 0, u) and (O, 0, Cl'), respectively. We have further supposed that 
the orientation of vorticity and microrotation are parallel. Under these 
assumption and in view of equations (1) and (2). The equations of motion can 
be written as 
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Y d 2 Q - 2n du - 4ni Q = 0 

dy 2 (jy 

On introducing the following non dimensional quantities 
’I = u , fi = Q 

h (Ah2/2 n) ~(A~h72jo 

X. = h 


(4) 


o 


a 




(5) 

( 6 ) 


C y (n + n,)/ 4n njM 
the equations (3) and (4) reduced to 
-1 + % (1+a) u" + a Q' = 0 
o'” - X £2' - X = 0 

where premes denote differentiation with respect to n. 

The equations (5) and (6) are solved under two limiting boundary, conditions. 

>) No spin boundary which corresponds to no-slip condition. Q (o) = n (1) 

0. This is based on the argument that fluid-solid interface with interaction 

IS so strong that the microstructure does not rotate relative to the plate. 

Weak limit of influence of microrotation on the surface given by Q’(0) = 

n'(l)=o. This is based on the fact that the corpuscles rotatate on the 
wall. 

In view of these two conditions and no-slip velocity at boundary, 
u(0)=u(l)=o, the solutions are. 

i2s - cosec h X sxnk X ti - n 

cosec h X [l-r|+Ti cos hX-cos hXTi]-Ti (l-n) (g) 


iO 


X(l + a) 


tivv - cosec h X [ cos h X T] - cos h (I-t]) x] 


(9) 


Uw 


2a [t)+(l-n)cot h X+{Ti-sln h Xti-cos h(l-Ti)X}cosec h X}]-rt(l-Ti) 

(10) 
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subscripts s and w denote fho . i 

■ ^ r . j ^ Strong and weak limits of 

mteracfon of flu.d-solid interface, respective^. 

Discussion : 

The nu^er.cal values of microrotation and velocity for different field 
parameters are obtained and tHeir variations in tHe flow regions are 

encountered from tHe taHie .e oHserve tHat ail tHe values oHtained are negative. 

This shows that backward flows occur wb: u • • ' 

obvious because the 

pressure gradient under whose influence tho fl 

. . te, be maintained, 

. , ' P'^“ f-- the high pressure 

to low pressure direction. From table I we infer that H • ii 

inter that having all its physical 

properties to the constant, when the width of the cha , ■ 

^ channel increases then the 
microrotation Q increases for strone limit ft 

limit ■ t . * ^ interaction while for the weak 

limit of interaction an opposite behaviour is seen Below tu 
w Below the central resrion 

From table II we see that fee. h 

iee mat tor small a velocitv ficM : 

u.i,„ ^ , velocity field increases with X below 

the central region and above it a reverse trend f 11 

of inter. ,' • ■ • 'whereas for weak limit 

whoh increase in velocity with ^ is recorded throughout the 

region. Also, when a increases velocity field decreases. 

9 

Variation of micrororational field 


n 

X.= l 

Qs 

A,= 3 

0.0 

-0.0000 

-0.0000 

0.1 

-0.0,14 7 

-0.0696 

0.3 

-0.0408 

-0.1975 

0.5 

-0.0565 

-0.2874 

0.8 

-0.0442 

-0.2543 

0.9 

-0.02 65 

-0.1607 

1.0 

-1-0.0000 

-0.0000 



X=5 

A.= l 

Qw 

A.= 3 

X=5 

-0.0000 

-0.4621 

-0.3017 

■0.1073 

-0.0929 

-0.4642 

-0.3138 

■0.2182 

-0.2713 

-0.4786 

-0.3908 

-0.3383 

-0.4184 

-0.5000 

-0.5000 

-0.5000 

-0.4322 

-0.5300 

-0.6545 

-0.7305 

-0.2935 

-0.5357 

-0.6861 

-0.7817 

-0.0000 

-0.5378 

-0.6982 

-0.8026 
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Variation of 

velocity 

field { a 

=0,2) 


r| 


Us 



Uw . 



x=i 

X=3 

X=5 

X=1 

X=3 

X=5 

0.0000 -0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

0.1 

-0.0760 

-0.0804 

-0.0835 

-0.0692 

-0.0811 

-0.0847 

0.3 

-0.1766 

-0.1846 

-0.1916 

-0.1609 

-0.1908 

-0.2006 

0.5 

-0.2091 

-0.2147 

-0.2219 

-0.1903 

-0.2294 

-0.2387 ; 

0.8 

-0.1324 

-0.1300 

-0.1312 

-0.1201 

-0.1445 

-0.1516 

0.9 

-0.0742 

-0.0712 

-0.0705 

-0.671 

-0.0808 

-0.0847 1 

1.0 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

i 

-0.0000 ji 
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Chapter Sr99 


St 3 day flow of blood in narrow tube through 
micropolar fluid model with slip at the wall. 


introduction : 

Blood, a class of fluid has its component as plasma, platelets, erythrocytes 
and leubocytes. But most of its rheological properties depend on RBC which 
in normal condition occupy some 40-45% fraction of total blood volume. The 
general evidence is that the flow processes do not obey the same rule everywhere 
as they vary from small to large diameter tube. In absence of our knowledge 
of cellular geometry and the flow inside the cell we assume, it a rigid particle 
of finite dimension suspended in base fluid plasma. In narrow tubes (20'500pm; 
1 H m = 10'^ m) we may propose the micropolar fluid theory as a model to 
explain the flow parameters characterizing the blood behaviour. Several authors 
like Klime (1968), Brunn (1975), Ariman (1974), Cowin (1972) Erdogen (1980) 
and Chaturani (1984) have proposed micropolar fluid theory a suitable 
candidate for the same. They have used different boundary conditions at the 
wall. 

Further, it is found to be more realistic to consider the shear viscosity 
of blood as different from that of solvent viscosity. In the present section we 
have studied the poiseuille flow of micropolar fluid as a model of blood flow 
in 40 ^m diameter tube and for 40% RBC concentration. 

/' '■ the boundary we have introduced a slip in the axial velocity and a 
partial rotation depending upon the wall effect parameter S and fluid vorticity. 
From the analysis we have observed the the introduction of slip reduces the 
apparent viscosity of the suspension and increases its axial velocity whereas 
it has no effect on particle rotation. Effect of S (0 < S ^ 1, a parameter 
depending on S and concentration) on different flow parameters are seen and 
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by comparing the results with the experimental value we have tried to estimate 
the value of S. 


j^athematical Analysis : 


Fully developed steady incompressible flow of micropolar fluid through 
a rigid circular tube of radius R has been considered^ As a model for the flow 
of blood in narrow vessels, Neglecting body forces and body couples, we write 
the equation of motion in the form ; 

(f + O 1 ^ a v ) + 2 X ^ Crw) - ^ = o CD 

dz 


3r 




r 2>r 


S' r 1 (rw) 


br -r ar 


. 2 X - A xw = 0 

at 


C2) 


where (0, 0, vCr)l is the velosity. wCr) the total angular velocity; r. z 
the radial and axial co-ordinates, p the pressure, u the shear viscosity, t the 
rotational viscosity and t the rotation parameter. 

Boundary conditions for the flow arc 
at r = 0; v and w remain finite 

at r = R; v = - a ^ w = ^ ; (0 ^ s ^ D 

^r 2 2r 

where a - the slip parameter, s - the wall parameter. 

Solution of equations (1) and (2) in view of (3) are 
V = K r 2 [l-y^ - 2n i ( lo CD - loC^y))! + 2 x Cl-^'S) 


C3) 


4 fi 


1 ii CD 


w = 1C R 


4 

where = 


y - i II C^y) 


C4) 

C5) 


ii CD 

y = r/R. ^ ' 


dz 


r (^ + t) 





h 


‘I 


■ r 
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H + T 


, i = 1 - s ; (0 < S < 1) 

y,^- 

1- T1 S 


is the modified Bessel function of order n. 

From equation (4), the flow rate Q, apparent viscosity Hg and velocity 
Va at tube axis are obtained as 

Q = 71 K [1- 4ti s { Iq (?.) - 2 } + 4 H/ (l-ri s)] (6) 

8 (.1 , ^ 

Ha = M [ 1-Vn i ( IqM - 2 ) + 4 VH , (l-T, "s)] ^ (7) 

X l^ix) 

Vg = KR2 [ 1- 2 ti i ( Io(^)-l ) + 2 V (1- T1 i)] (8) 

”47 ^ 



Discussion : 

From equation (4) and (5) we obserc that slip at the boundry affeels 
the axial velocity whereas the particle rotatoin is unaffeted. 

The boundry condition parameter of Brumn (1975) analysis is related 
with S by the relation B^ = Jj ^ 0 for the discussion we have considered 

two case. In case (I)- we suppo^ n = Hs + t and in case II we suppose 
H = Hj where Hs the solvent viscosity. We see that the results obtained for 
case I are in conformity with the experimental value, wheras the results of 
case II are in parity to the results of Erdogen (1980) analysis which on 
comparing with the experimental value give unsatifactory reading for 40% 
concentration and 40 pm tube we have used. 

Hj = 1.2 Cp, X = 0.98 Cp, r = 12 x lO'® gm cm/sec. from Ariman (1974) 

analysis. 

For exerimental verification we have used the data of Bugliarello and 
Sevilla (1970). 

For different values s, vp. the theretical values of Ha> ''a- 

cTZ. 

w(r) are shown in tables I - V. 





For v|< = 0.02, the apparent viscosity varies from 2.02 to 2.8 where 
as for V); = 0 this range is from 2.18 to 3.109 which is greater than the pervious 
one. From the analysis of Bugliarello and Sevilla, the cxperimantal value of 
for 40% concentration could be read any where between 2.1 and 2.5. To 
discuss the other parameters we have chosen the value of |j,a arbitrarily equal 
to 2.5 Cp. We are not known with any precise experinental value of i. It can 
be estimated only from equation (7) as in it other quantities except i are 
experimentally known. The values of s lying between 0.35 and 0.66 give better 
theoritical result. We can approximate s = 0.5 because for this. 

Value of s there are small errors in Hg and Vg when compared with 
experimantal values. 

Variation of pressure gradient and axial velocity with respect to \f and 
s are reported in table II. 

Variations of velocity and rotational. Profiles for different s are shown 
in tables III, IV From tables we observe that when s increases, particles axial 
and rotational veloities decrease. If s is assumed to be the measure of 
concentration, then it is seen that for dilute suspension, particles rotation 
increases with tube radius (a case of Newtonian fluid) . But for sufficently large 
value of s (s > 0.85) the rotation at first increases, attains its maximum value 
and then decreases near the wall. For maximum concetration (i = 1) the particle 
rotaion is maximum of r = 0.75 which has been experimentally observed. 

l^eferences : 

Ariman T, Turk, M.A. and Syluester N.D., Steady and pulsatile flow of blood. 
Jour. Appl. Mech. ASME, 41, 1, 1974. 

Burmn, P., The velocity slip of polar fliuds. Rheologica acta. 14, 1039, 1975, 
Bugloarello. G, and Sevilla, J, velocity distribution and other characteristics of 
study and pulsatile blood flow in fine glass tubes. Biortheology, 7, 85, 
1970. Chaturani P and Biswas D. A comperative study of poisewille flow 
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fSovA acte. 23. 435, 1984. 

Cowin, S.C., on the polar fluid as a model for blood flow in tubes. Biorhcology, 
9. 23, 1972. 

Erdogen, M.E. polar effects in apparent viscosity of a suspension. Rheol Acte. 
9. 434, 1980 

Kline, K.A. Allen, S.J. and Desilva, C.N., A continuum opproach to blood flow. 
Biorheology, 5, 111, 1968. 


^able 9 

Variation of and for different values of S 


1 



s 


0.591 

0.591 


<€able 99 

V and Vg with respect to y and s 


-dp 


65.31 

64.27 


0.048 

0.097 
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stable V 

jriation of appar ent viscosity and axial velocity 
for different values of S at constant xy = 0.02 


30 

50 

00 

,66 

)00 


1.11 

1.31 

1.42 

1.56 



5.66 

5.35 

4.44 

4.03 

3.22 


0.217 

0.183 

0.152 

0.152 

0.119 



study of dispersion processes in blood flow in 

narrow vessels 


In biological sciences, the study of diffusivity of neutricnts, metabolic 
products, drugs and other solutes is of utmost importance. Specially, many life 
giving materials mixed in the blood reach to the different parts of the body 
by the process of diffusion. Taylor (1953, 54) studied the dispersion process 
in Newtonian flow and discussed the effective dispersion coefficient with respect 
to the average speed of the flow, the radius of the tube and molecular diffusion 
coefficient. Many authors as Patel and Sirs (1983), Federspial and Fredberg 
(1988), Rudraiah et. al (1986) have studied dispersion processes by taking 
different flow models. Shukla and Gupta (1981) studied the Taylor dispersion 
in Bingham Plastic fluid surrounded by Newtonian peripheral layer, Scottblair 

Charm anH ICu rland (1968). Bueliarcllo and Sevilla (1970) 
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j^athematical analysis : 

Consider an incompressible steady viscous flow of casson fluid through 
a circular tube of radius R. During motion, the viscosity and yield stress of 
fluid vary along the radial direction. The constitutive equation in one 
dimensional form of casson fluid is ; 

K’^"(r) + T] "^(r) x > ly K(r) 

> 

Y = 0; T < ty K(r) CD 

where nC?") 's the viscosity, XyK(r) is the yield stress, y 's the strain 
rate. r|(r) and !<(r) are assumed to be decreasing function of r. The region 
0 < r < R is divided into a central region 0 < r < of casson fluid with a 
plug radius Rq. 0 < r < Rq < and by a peripheral region of another casson 
fluid (Rj < r < R). The function Ti(r) and K(r) are assumed to be 

r| (r) = 111 ! KCr) = 1 when 0 < r < Rj 

XI r) = 11 2 i ^ ^ *^Cr) = K < 1; when R^ < r < R (2) 

Where and ri2 viscosities of the central and the peripheral 

fluids. In case of blood, the peripheral layer is a Newtonian fluid and thus the 

function K(r) = k = 0 in R^ < r < R. 

For one dimensional steady laminar flow in cylindrical co-ordinate system 


a 

ry 
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In view of equations (4), (5) and boundary condition at r-R, the 

integration of equation (1) gives 


/ -r 2 dp W r'y dy ^ yo r 

K(y) dy 

2yo'‘ C 

K«(y) y« dyl 

dx >' Ly lf(y) K(yO) y 

TiCy) 

K«(y„) y 

hCy) J 


: yo < y < 1 (6) 

Where y = r; yo = Ro> 

R R m 

For y = y^, the plug velocity is calculated from equation (6) as 

-R2 dp^r^'ydy yo K(y)dy 2yo^ K^(y)y^ dy; 

Wrii dx' L yo iiCy) K(yo) yo n(y) K’^Cyo) yo il(y) 

0 < y < yo 

The average velocity V, of the fluid is obtained as 

V = 1 f 2nr Vx dr = - y^ dVx dy 

nR2 o ° dy 

Usi.ig equations (6) and (8) in above equation, we get 

V = 1 ^ R2 dp\ r |'y3dy yo j^'y^ K(y) dy K'^Cy) dy 

2tii dxHyo Ti(y) K(yo) yo nCy) 


Now, the equation for concentration is 
^c + dc = |r D(r) 2) c 1 


K’^Cyo) yo 


(9) 

( 10 ) 


Where D(r) is moleculer diffusion coefficient assumed to vary symmetrically 
along r- d irection. 

Now introducing the non-dimensional quantities 

0 = t, t = L. ^ = X - V t. y = r. D - D 


Where L and D are the typical values of The length along the tube and 
molecular diffusion coefficient, respectively. Equation (10) is transformed to 
give flow of solute relative to a plane moving with the mean speed of the flow 




m 

I ’“i- 


lac + (vx ■ V) a c « D 1 ^ iy 0 he 


t 2te 


L R2y 2»y 


( 11 ) 


Applying Taylors limiting conditions ( ^ is independent of y and _^*o), 


equation (11) gives 


y ay 

where 


D y he' -= F g (y) 


(12) 


ydy yo f K(y) dy ly^ 


y^ K’^(y) dy - G; 


y'o ^(y) K(yo) yo n(y) K’^Cyo) yo ^(y) 

; 0 < y < yo 

g(y) = f ydy yo f K(y) dy lyo"^’ r y"^ K’^(y) dy - G; 


yo n(y) K(yo) y^ n(y) K’^(yo) yo ^^y) 


yo ^ y ^ 1 


^^y3dy y^ T ^ y^K(y) dy Zyo"^ f y^'^ K’^^y) dy (14) 


^0 "^y) K(yo) yO “-n(y) 
/- 1 dp. \ a c 


K’^(yo) yo ^(y) 


(15) 


2tii d^/ DiL2 3^ 


we obtain 


Int;,_<,rating equation (12) and using conditions c=0 at p=0; ac=0 at p = l. 


C=Co - F M (y) 
where M(y) = 


Dy y 


y g(y) dy> dy 


(16) 

(17) 


The average volumetric flow rate ^ is given by 
_ I _ 

Q = 2 Jy (V^ -■ V) c dy 

Using equations (6), (8), (9) and (16) in the above equation we get 

Q = — ,2V2 r 2 fy M(y) g(y) dy 

' 

□ iL r r'y3 dy y^ f' y^Kiy) dy - 2 yo^^ K^Cy)dy 

2'nCy) ”^K(yo) yo K^Cyo) yo ”^(y) 

(18) 



37 


Compare equation (18) with the Picks law of diffusion 
-D* 

d X 

The solute disperses relative to a plane moving with the menan speed 


of flow with an effective dispersion coefficient D*, given by 


D* = 2R2v 


Ja 


o ■ 


y M(y) g(y) dy 


fr' y’dy 

+ yo 

Jy2K(y) 

dy - 2yo’^ 

jr'y5/2K%(y) dy' 

[yi V(y) 

K(yo) 

yo ”n^y) 

K^(yo) 

yo ~^y) 


(19) 

The expression for D* is applicable for any general functions 
^(y). D(y) and K(y). 

Now, for the effect of a peripheral layer around a casson fluid, the 
boundary condition for the functions Ti(y), D(y) and K(y) may be assuemed 
as follows : 

Ti(y) = fl i 0 ^ y < yi = 

R 

n2 = 1 yi ^ y ^ 1 

^1] 

D(y) = [ 1 ; 0 s y < yi 

^2 = J yi - y - ^ 

Pi 

K(y) = Jl ;0<y<yi 

\o ; yi s y < 1 

where yl, and Dj^ are the radius, viscosity and .molecular 

diffusion coefficient of casson fluid respectively and ^2* *^2 peripheral layer 
fluid. 

Using equations (20), (21) and (22) in equations (13), (14), (17) and 

(19) we obtain 

D* = r2 v2 H 


( 21 ) 

( 22 ) 


D 


a 
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Where 

H * Sq (S]^ - S2/v) 

[ yi' • yo‘^11 {yo""(‘ - '^yl’ - 12 

21 

Sq * 1 - (I'h) C2-yi^) yi^ - yo^^^ ■*■ ^yo>'iH 

21 3 

- ISyo"^’ (7-3yi2) 

21 a 

Si = yi^ ioge yi ■ ^ - yj^ ^ 

2 4 8 

$2 * Si + 3 , 

8 

Discussion : 

For fixed mean speed of the flow, the effects of yo and v, namely the 
viscosity, yield stress and molecular diffusion coefficient on H are seen from 
the tables (I), (11) and (III). We find that the magnitude of H decreases as 
the yield stress and molecular diffusion coefficient increases but increases with 
respect to the coefficient of viscosity. 


(cable 9 

Effects on y., v and y^ on H, ( v = 1.1, yj - 0.90) 
H 


yo 

p =0.20 

0.40 

0.60 

0.80 


0.30 

-0.0480 

0.0873 

-01316 

-0.1710 


0.05 

-0.0321 

0.0680 

-0.1034 

-0.1360 


0.09 

-0.0219 

-0.0450 

-0.0689 

-0.0928 
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^able S7S7 




(v 

= 1.3, yj = 0.90) 


H 



yo 

H=0.20 

0.40 

0.60 0.80 

0.03 

-0.0401 

-0.0857 

-0.1292 -0.1679 

0.05 

-0.0315 

-0.0668 

-0.1015 -0.1336 

0.09 

-0.0215 

-0.0442 

-0.0677 -0.0904 


^ab/e 999 

v = 1.5, y^ = 0.90 


H 


0 

yo 

H = 0.20 

0.40 0.60 

' 0.80 

0.03 

-0.0395 

-0.0846 -0.1275 

■0.1657 

0.05 

-0.0311 

-0.0659 -0.1002 

•0.1318 

0.09 

-0.0212 

-0.0436 -0.0668 

■0.0892 


Hence inview of equation (23) and above discussion it can be seen that 
for a given mean speed of the flow, the effective dispersion coefficeient D* 
decreases w'ith increasing and v while it increases with respect to p,. 
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Chapter V 


A theoretical Analysis to the 
Effect of Wall Cayer thickness on l^lood Kheology 

introduction : 

oka (1967) and Jones (1966) have explained the Coplay -Scott Blair result of 
the loss of apparent blood viscosity in a glass tube when the wall surface is 
coated with fibrin solution by introducing a slip velocity of the surface. 
Whitemore (1967) and Hyman (1973) have suspected the view and opine that 
slip at the solid surface has no evidence in practice and therefore is not a 
viable alternative to explain this reheological behaviour of blood. From the 
several tudies it has been established that for blood flowing in narrow tubes 
(500 ixm to 50 iim), a cell free plasma zone exists near the wall and the extent 
of this depends on many rheological parameters. From the analysis presented 
here, we see that apparent viscosity of blood decreases as the thickness of plasma 
layer increases. But there is no precise quantitative evidence that why this layer 
increases when its molecules come in contact with fibrinised surface of the wall. 
In one of his result Oba (1981) has contended that the thickness increases 
due to electric change on fibrin molecules. 

I Mathematical Analysis 

Two fluids model of blood with a central core region of radius R and a peripheral 
layer of thickness 5 = R - RO, both satisfying casson equation 

tVi = Xy'^' + 


has been considered, where t represents shear stress, q the casson 
viscosity, y the strain rate, Xy the yield stress and R the tube radius. 
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The equation of motion for steady incompressible laminar flow 
^ p = 1 5 (r 

d z r br 

When integrating with respect to the radius co-ordination r, and assuming 
that the stress at the axis be finite, gives 

X = d p r C3) 

2 


d Z 


Where p is the pressure, z the axial co-ordinate, x^z the shear stress normal 
to r in z direction. Substituting equation (3) in (1) and assuming that 


'^rzCZ) 


= t 


rzCl)» ^(1) “ ^^(2) at r = Rq 


and Vq) = 0 at r = R 

the velocity in peripheral and core regions are obtaineci^ as 
= rxr [1- r2 - Pci) 

3 V R5/2/ 


2q 


1 L 


R2 


V(2) - 


2ti 


^1- r2 8 x3/2 /I- r3/2Y 2 P 2 X fl* r 

3 ' 




R 3^: 
^0 - 


(4) 

(5) 


+ Rxr [( 1 -x 2 ) P^)’^ (1-x3/2) + 2pi (1-x)] 

3 


2tii 


Where, x = 1 - 6 . Pi = Xy^), P 2 *= x„o, xp is the shear stress at the 


Xy2> '■R 
Xc 


R Xr IR 

wall and suffices (1) and (2) denote the values of the corresponding parameters 
in peripheral and core regions respectively. 

Plug flow exists wherever the shear stress does not exceed yield stress 
and can be obtained by putting r = P 2 R in equation (5). From (5) we may 
derive Jones result of one fluid model of blood flow. 

We calculate apparent viscosity xia in the form as 


1 = 1 1 

: 1- 16 P2’*"* -»• 4 p 2 - 1 Pz^^i 

+ 46 

;i+Pi - 2Pi’^ 

- I+IS2 

Ti2 

[ T~ 3 21 

R 

L m 

nz 





(6) 


i 


a 

y 
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Wc observe that the loss in is caused due to increase in 5 
(as Ti 2 > 1 + P 2 ■ 
ni 


1 + Pi - ZPi”^ 


Results are discussed in the form of following two cases 
Case I. When p^ = 0 and the relation tij = tjj (1 - a (j)) holds, where 
^ is hematocrit function and a is experimental constant. Then 

1 = 1 (1 - a (t.) [1 - 16 + 4 P 2 p 2 ^] C7) 

’ll 7 3 21 

+ 1 [(1 - a<lj) (2P2’^ ' P2^ + a <j)] 45 

’ll 


R 


The above equation (7) shows that for any desired concentration (|i, tig 

decreases with 5. 

Case II. When Pj = 0 = P 2 

In this case the above theory gives the result for Newtonian model in 
which the fluids of different viscosities r\-y and n2 both regions but have 

O' ^ 

the same constitutive stress-strain relation of Newtonian structure. The appar- 
ent viscosity Tig is calculated as 
1 = 1 [1 - (1- 45) a 

ila 'll ^ 

Equation (8) shows that for 5 = R/4. Ha = ’ll and when 5 > R/4 then 
tig decreases with increasing 5. 

This supports the copley second result in which the experiment has been 
performed for plasma fluid alone. 
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theoretical Model of Steady l^lood Jlow 
in flat row Vessel 

introduction : 

Mathematical modelling of blood flow in different vessels under various 
physiological conditions have drawn the attention of several workers. Since 
poiseuilles empirical euqation of finding out the apparent viscosity of viscous 
fluid in fine glass capillaries, many workers attempted to find the apparent 
viscosity invivo and invitro, which is an important parameter affecting the blood 
rheology. Fahracus and Lindquist found that apparent viscosity of blood 
decreases with the decrease of tube diameter below 500 ^m down to 7 ^im. 
Recently, Kiani and Hudetz (1991) developed the model assuming no yield stress 
behaviour of blood in these vessels. It has been recognised by several 
experimental observations and viscometric datas that blood possesses consid- 
erable amount of yield stress in narrow vessels. Scott Blair (1959), Charm and 
Kurland (1968), bugliarello and Sevilla (1970), Oba (1981, 89) have proposed 
casson constitutive equation to analyse the blood behaviour in narrow vessels 
at low and high shear rates. In the present analysis we have assumed that the 
fluid in core region satisfy casson equation and in marginal layer region satisfy 
Newtonian equation. Apparent viscosity of blood has been determined as a 
function uf yield stress, vessel diameter and peripheral layer (PPL) thickness. 
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jdathamatical /Inalysis : 

Incompressible fully developed laminar flow of blood is considered in a 
circular tube of small diameter. Two-layer flow model is assumed. Core fluid 
follow the casson equation. 

(1) 

and plasma peripheral layer follows the Newtonian equation 


t t^p 


( 2 ) 


Where x is shear stress, Xy the yield stress, y the shear rates, casson 
core viscosity, i-ip Newtonian plasma viscosity. 

Steady incompressible laminar flow equation 

ap = 1 (r x,.p 

dZ r dr 

after integration yields 


Tz 


= P*' 
2L 


(4) 


Where 

z, r = axial and radial measurements 
t|,j, = shear stress normal to r in z direction 
, p = pressure drop along length L of the tube. 

i Al the interfacial boundary between core and marginal fluids, shear 

stresses are equal and at the wall surface noslip condition are assumed. 
Substituting (1) and (2) in (4) we get 

Vp = P (Rw2 - r2); < r < Rw (5) 

4L pp 

= P (R<-^ - r^) - 4 /PTy^ 1 

4Lp,. 3 V Pc 

^ + Xy (Rc - r) + P (R^2 . r^ 2 ) ; 0 ^ r £ Rc (6) 

Where Vp denotes plasma velocity and Vc denotes core velocity. 

Plug velocity is obtained by putting r*Rc in equation (6) as 
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V, - P (Rc^ - R,2) + V (R, ^ R,) 

4Lnj. 3\x^ V 2L/ 

+ P (R^2 . r^2) 


(7) 


4Ln, 


The total flow rate Q = Q + Q^. + Q| 

g R 

or Q = f 2:irdr . Vp + J Zurdr . V^. + 7t R|^ . V| 


R, 


R 


or Q “ PkR 


w 




+ 4 / 


8 V R. 




R 


w 


R. 


'W 


Wh-re plug radius R| = 2 Lt 


( 8 ) 


The apparent blood viscosity Pgpp is obtained as 

a - 25\^ - 16 /I- 25 , <jV4 



I 


- 


d^ 


(9) 


a = R|, 5 = R^ - R^, d = 2R, 


R. 


lipto first order smallness of 8, equation (9) reduces to 


Papp 


p.p »F(a) + 85 (1 - ^ip <j)(a)}l 


( 10 ) 


Pc d \ p, / i 

Where F(a) = 1-16 + 4 a - 1 

7 3 21 

(l)(a) = 1-2 4- a 


a 
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As a is a measurement of yield stress, therefore at high shear rates 
it takes small value but not zero as chosen by Kiani and Hudetz (1991). For 
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constant F(a) and <t»(ot), apparent viscosity decreases with the increase in 5. 
Fahracus Lindquist effect can be observed from equation (10). Charm and 
Xurland have shown that 8 increases with decreasing a. Wc have taken the 
datas of Charm and Kurland and shown the variation of with a at constant 

•"app 

in table I. Assuming that the casson viscosity is related by general equation 


(1-KH) -1 (11) 

Where H is a function of hematocrit and K is an experimental constant, 
we can obtain ^gpp at different concentration of hematocrit. 


^A^£E 9 

[charm and Kurland (1968), Hp = .0125 poise. 


= .0259 pose, d = .01, H = .448] 


^app 

a 

F(a) 

f(a) 

2.7790 

.0005 

.9495 

.9559 

2.8293 

.0010 

.9291 

.9378 

3.0584 

.0050 

.8451 

.8636 


a 
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